Indirect exchange interaction between magnetic adatoms in a monolayer M0S2 
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We study the Ruderman-Kittle-Kasuya-Yosida (RKKY) interaction in a monolayer M0S2. We 
show that the rotation of the itinerant electron spin due to the spin-orbit coupling causes a twisted 
interaction between two magnetic adatoms which consists of different RKKY coupling terms, the 
Heisenberg, Dzyaloshinsky-Moriya and Ising interactions. We find that the interaction terms are 
very sensitive to the Fermi energy values and change dramatically from doped to undoped systems. 
A finite doping causes that all parts of the interaction oscillate with the distance of two magnetic 
impurities, R and the interaction behaves like R~ 2 for the long distance between two localized spins. 
We explore a beating pattern of oscillations of the RKKY interaction which occurs for the doped 
system. 
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I. INTRODUCTION 

Two dimensional (2D) materials can be mostly exfo- 
liated into individual thin layers from stacks of strongly 
bonded layers with weak interlayer interaction. A fa- 
mous example is graphene^ and its analogs such as boron 
nitridfP. The 2D exfoliate versions of transition metal 
dichalcogenides 3 exhibit properties that are complemen- 
tary to distinct from those in graphene. M0S2 ( the 
polytype is 2H-M0S2 has trigonal prismatic coordina- 
tion) is a hexagonal crystal layered structure with a co- 
valently bonded S-Mo-S hexagonal quasi-two dimensional 
network which does not have inversion symmetry, packed 
by weak van der Waals interactions. The monolayer of 
M0S2 has provided a new material with a peculiar struc- 
ture for the charge and the spin interactions^. There is a 
transition from an indirect band gap of 1.3 eV in a bulk 
structure to a direct band gap of 1.8 eV in the monolayer 
structure^. This direct band gap opens the possibility of 
many optoelectronic applications^. The electronic struc- 
ture of M0S2 also enables valley polarizatiorPED since 
both the conduction and valence band edges have two 
energy-degenerate valleys at the corners of the first Bril- 
louin zone and valley selective circular dichroism arising 
from its symmetry. 

Ruderman-Kittle-Kasuya-Yosida (RKKY) interaction 
is an indirect interaction between nuclei of transition 
metals or magnetic impurities, mediated by the conduc- 
tion electron s ! 10 ! 12 ! 13 !. This interaction is directly related 
to the spin susceptibility of the host metal. Rudermann 
and Kittel first suggested^ that the spin oscillatory in- 
teraction in metals could provide a long-range interaction 
between nuclear spins and explored the effective magnetic 
interaction between nuclei of transition metals and thus 
described the reason for broadening of spin resonance line 
diagram. Due to the fact that the RKKY interaction is 
originated by the exchange coupling between the impu- 
rity moments and the spin of itinerant electrons in the 
bulk of the system, a spin splitting owing to spin-orbit 
coupling, must have opposite sign at the two valleys by 
the requirement of time-reversal symmetry, is expected 



to influence directly this interaction. Since the two in- 
equivalent valleys in the monolayer M0S2 are separated 
in the Brillouin zone by a large momentum, in the case of 
the absence of short-range interactions, intervalley scat- 
teringpS should be negligible and thus the valley index 
becomes a new quantum number. Therefore, manipu- 
lating valley quantum number can produce new physical 
effects. 

In the monolayer M0S2, the highest energy valence 
bands and the lowest energy conduction bands, at the 
K and K' valleys in momentum space, are mainly of 
molybdenum d— orbital character and spin-orbit interac- 
tions splilJSIthe valence bands by 2A ~ 150 meV leading 
to a spin polarization of the valence band. A minimal ef- 
fective band model Hamiltonian of M0S2 is found by Xiao 
et aP^ with parameters based on first principle calcula- 
tions. To first order in momenta, the first part of Hamil- 
tonian describes the dynamics of massive Dirac fermions 
which are studied from the graphene committees. It is 
commonly known that the RKKY interaction in massless 
graphene is quite differentP^HIZI from that of a shcrodinger 
2D electron liquid. The second part of the Hamiltonian 
illustrates the spin-orbit interaction which leads to cou- 
pled spin and valley physics in the monolayer M0S2. The 
RKKY interaction usually yields a Heisenberg coupling 
refereing to a parallel or antiparallel coupling of local- 
ized spins. If spin of conduction electrons processes, it 
can be possibly to produce a noncollinear Dzyaloshinsky- 
Moriaya ( DM) coupling of localized spins. We thus use 
this model Hamiltonian to calculate the exchange inter- 
action between localized spins on the M0S2 and show 
that, in particular, combination of the spin-valley with a 
massive Dirac-like spectrum can mediate a much richer 
collective behavior of magnetic adatoms. 

The RKKY interaction in 2D electron gas in the pres- 
ence of a Rashba spin-orbit coupling has been studied 
by Imamura et alP^ and they showed that rotation of 
the spin of conduction electrons causes a twisted RKKY 
interaction which consists of three different terms, the 
Heisenberg, DM and Ising interactions. Moreover, the 
RKKY interaction between localized magnetic moments 
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in a disordered 2D electron gas with both Rashba and 
Dresselhaus spin-orbit couplings has been studied^. The 
disorder-averaged susceptibility leads to a twisted ex- 
change interaction suppressed exponentially with dis- 
tance, whereas the second-order correlations, which de- 
termine the fluctuations of the RKKY energy, decay with 
the same power law as in the clean case. 

In this work, we calculate the RKKY interaction me- 
diated by spin-orbit coupled massive Dirac fermions in a 
monolayer M0S2 using the Green's function method. We 
show that the interaction consists of three terms with- 
out having the Rashba or Dresselhaus interaction and 
show particularly that the interaction behaves like R~ 2 
for the long distance between two localized spins. The 
DM coupling leads to a rotation of the spinors of mag- 
netic moments with respect to each other. In addition 
a beating pattern for the interaction, in the cases where 
the monolayer M0S2 system is doped, is obtained. As 
we mentioned before, the model Hamiltonian can be re- 
duced to a massive graphene by discarding the spin-orbit 
coupling. Besides i?~ 3 dependence of the interaction for 
undoped massless graphene, we find that the interaction 
decays like Br 2 for a massive and doped graphene. How- 
ever, for a massive and undoped graphene, the interaction 
decays rapidly in a good agreement with those results ob- 
tained in Ref. [2"U] . 

The paper is organized as follows. In Sec. [II] we intro- 
duce the formalism that will be used in calculating the 
RKKY interaction from the Green's function. In Sec. IIIII 
we present our analytic and numeric results for the cou- 
pling strengths of the RKKY interaction in both undoped 
and doped monolayer M0S2 sheets. Section IV contains 



discussions and a brief summary of our main results. 



II. METHOD AND THEORY 

The Hamiltonian for a single layer M0S2 including two 
magnetic impurities can be written as % — Ho + Hint, 
where Ho is the noninteracting Hamiltonian for the single 
layer M0S2 and Hint is the interaction term between im- 
purities and spin of conduction electrons. In the mono- 
layer M0S2, the highest energy valence bands and the 
lowest energy conduction bands, at the K and K' val- 
leys in momentum space, are mainly of molybdenum 
d— orbital character and spin-orbit interactions splrP^ 
the valence bands by 2A ~ 150 meV leading to a spin 
polarization of the valence band. Fig. [T] shows the band 
dispersion of monolayer M0S2 in which spin splitting in 
conduction band is negligible however a spin splitting in 
the valance band leads to a spin shift. 

The Hamiltonian for M0S2 in a continuum model at 
K and K' points is written aJ^ 



Hq = at(rk x a x + k y a y) 
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where the index r is ±1 at valley K or K' , A is the direct 
band gap, A is the spin-orbit coupling constant, a is the 




FIG. 1: (Color online) Schematic picture of the dispersion 
relation in momentum space at valley K. The valance and 
conduction bands are separated by a large direct band gap, 
A and there is a large spin-orbit coupling, A, leading to a 
spin polarization of the valence band, spin up t an d down 
spin, 4. The dispersion relation is given by E(k) = \ts z ± 

VV fe ) 2 + ( A /2~ \rs z /2) 2 . 



lattice parameter and t is the hoping integral, tr's are the 
Pauli matrices written in pseudospinor ijy ~ ($^,$^) 
where c and v denote the conduction and valence bands, 
respectively and finally s z is the Pauli matrix for the z- 
componcnt of spin. 

Particulary, the noninteracting Hamiltonian in the 

spinor basis, = ( c l -p v \ *■> c k l> v k 1) can be written 
as 
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where 9' = tan -1 (k y / 'k x ) . The system in question here 
incorporates two localized magnetic moments whose in- 
teraction is mediated through a large spin-orbit coupled 
electron gas. The contact interaction between the spin 
of itinerant electrons and two magnetic impurities with 
magnetic moments Ii and I2, located respectively at Ri 
and R.2, is given by 



J, 



E 

i=i,2 



Ij • s(Rj 



(3) 



where J c is the coupling constant between conduction 
electrons and impurity, ( we set h = 1 from now on), 
s(r) = i 5(ri — r)<Ti is the spin density operator with 
Yi and Oi being the position and vector of spin operators 
of ith electron. 

The RKKY interaction which arises from the quan- 
tum e ffects is o btained by using a second order pertur- 
batioi jioii2i!3ii8| w j lic j 1 reads as 



W.RKKY = J c / J II Xii(R->R') H 



(4) 



where XylRR-') is the spin susceptibility in real space 
given by Xy (R, R) = E a6M X?/*(R, R'). ™ ^ h 
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X<f(R,R) - -^Tr[ [** deSmfaG*' (R,R',e) 
J 2?r J_ 00 

x a j G Pa (R',R 1 s)]] (5) 

Here a and /3 denote pseudospin degree of freedom, ef 
is the Fermi energy, trace is taken over the spin degree 
of freedom and a is the Pauli spin matrix. G a/3 (R, R'; e) 



is a 2 x 2 matrix of the single particle retarded Green's 
functions in spin space. In order to calculate the in- 
teraction Hamiltonian of Eq. Q, the form of the elec- 
tronic single particle Green's function, G(R, R';e) =< 
R|(e + i0 + - "H ) _1 |R' > is needed. To calculate the 
retarded Green's function in real space, its Fourier com- 
ponents in momentum space might be first obtained. 
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Generally, the energy dispersion of the noninteract- 
ing system can be obtained by the pole of the 
Green's function and it reads as E(k) = \ts z ± 

\J{atk) 2 + (A/2 - Xts z /2) 2 . The retarded Green's 
function are diagonal in the spin space, G cc (q,e) = 

t ) and G vv (q,e) = I ) 

V U P+-i 2 / V /3+-9 2 / 

where a± = (e + A/2 ± A)/(ai) 2 , a = (s - A/2) /(at) 2 
and p± = (e - A/2)(e + A/2 ± \)/(at) 2 . Taking the 
Fourier transform of the matrices, we thus find the real 
space Green's function matrix 

G CC (R, 0, e) = J d 2 q[e lK R e^ R G cc ( q (K)e) + 

e lK '- R e^ R G cc (q(K'),e)] (7) 

where q(K) or q(K') is the wavevector near the Dirac 
point K or K' , respectively. The integral can be simply 
calculated and it thus reads as 

G cc (R0,e) = 

/ e *K. V _ +e ,K<.R 5c+ x 

^ o e * KR ffc+ + e ,;K '' R g c J 1 j 

where 5 C ±(|R|) = -a ± K (y/-^± |R|)/(2tt) in which 
K o(x) is the modified Besscl function of the second 
kind. We can follow conveniently the same procedure 
discussed above to obtain G° c (0, R, e) = G° c (— R, 0,e). 
Moreover, the retarded Green's function G°„(R, 0,e) is 
obtained the same as the channel (cc) with replacing 
g v ±(\R\) = -aK (^^\R\)/(2ir) instead of 5c± (|R|). 

By inserting the retarded Green's functions in Eq. Q , 
and taking the trace over spin degree of matrices, the 
RKKY Hamiltonian simplifies to 



TiRKKY — JcXxxillxhx + Ilyhy) 

+ JciXxyih X h)z + Xzzhzhz] (9) 

which indicates that the RKKY interaction consists of 
three different terms: the Heisenberg like Ii ■ I 2 , Ising 
Iizhz and Dzyaloshinsky-Moriya (ii x l2) z terms. Ac- 
cordingly, the RKKY Hamiltonian can be written as 

7-Lrkky = JhHh + JiHj + JdmHdm (10) 
where Ji = P de^sm r/i and 

m = 4.9a+.ga-+2( 5 2 ++5 2_)cos((K-K')-R) 

a£ (c,v) 

m = J2 2(. 9a+ -.g a -) 2 (l-cos((K-K')-R)) 

a£ (c,v) 

VDM = £ - 2 (5a+-3a-) S M(K-K')-R) (11) 

a£ (c,v) 

The resulting RKKY interaction consists of three quite 
different interactions. The Heisenberg and Ising interac- 
tions favor a collinear alignment of localized spins. On 
the contrary, the DM coupling favors a noncollinear align- 
ment of localized spins. The same form of Hamiltonian 
has been obtained in 2D electron gas in the presence of 
Rashba spin-orbit coupling^ and in the topological insu- 
lators, on which magnetic impurities exhibit a frustrated 
RKKY interaction with two possible phases of ordered 
ferromagnetic phase and a disordered spin glass phased 
Importantly, the existence of Jdm in those systems is a 
consequence of spin changing term however in the mono- 
layer M0S2 system, the DM coupling is a consequence of 
the spin-orbit coupling and the existence of two inequiv- 
alent valleys separated by (K — K'). Consequently, the 
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DM term results in noncollinearity of two spinors inter- 
acting to each other. In addition, all coupling terms in 
the Hamiltonian for the monolayer M0S2 system depend 
on the direction between magnetic impurities. 



III. NUMERICAL RESULTS 

In this section, we present in the following our main 
results for the numerical evaluation of the RKKY ex- 
change coupling in massive Dirac fermions by analyz- 
ing the above calculated integrals of Jh, J dm and Jj. 
We sePI a = 3.193A, A = 1.66eV, 2 A = 150meV and 
t = l.leV in all our results presented in this Section, we 
also consider two valleys at K — (27r)(l, a/3) / (3\/3a) and 
K> = (27r)(-l,x/3)/(3V3a). 

At finite value of the Fermi energy, electron or hole 
doped, more complicated behavior of the RKKY cou- 
pling can be occurred. In this case, the behavior of all 
Ji interactions are determined by a superposition of two 
modified Bessel functions with two different periods de- 
termined by arguments in f3±. As the result, we observe 
that for doped M0S2, oscillations of Ji exhibit a beating 
pattern. Fig. [2^a) shows this beating behavior of inte- 
gral, Jr scaled by J,?/(27r) 3 a 2 as a function of the im- 
purities distance along x direction ( zigzag direction) for 
the Fermi energy which intersects the conduction band 
at £p = A/2 + A. Fig. [2jb) and (c) show the similar 
behavior for integral J dm and Ji, respectively as a dis- 
tance along x direction. Our numerical results show that 
all coupling interactions decay as R~ 2 at finite Fermi 



energy. As it is obvious from Eq. (10), in a long-range 



regime, Ji > Jh > J dm and moreover the coupling J dm 
and Jj vanish when A — > 0. Therefore, the existence of 
DM and Ising interactions is a consequence of spin-orbit 
coupling without having Rashba or Dresselhaus interac- 
tions. Besides, since g c ± are coupled to the Dirac points 
differently as diagonal elements in Eq. ([8]), the DM cou- 
pling is nonzero as long as one has the spin-orbit coupling 
together with two inequivalent valleys. It should be no- 
ticed that there are some locations along x direction for 
which J dm 3> Jh, since they have two different periods 
and this leads to a strong rotation even at long distance. 

Notice that in the M0S2 system, the RKKY interaction 
is direction dependant which is in contrast to the results 
obtained in the topological insulators or 2D electron gas 
systems. Since the RKKY interaction is a weakly inter- 
action and decays as R~ 2 in the M0S2 system, it would 
be good to look at its short-range behavior. Fig.[3j shows 
the angel dependance of Ji in contour plot scheme for a 
short distance and for ep = A/2 + A. Here, the sys- 
tem is electron doped and we assume that one of the 
spins is fixed at origin and another can be located on 
lattice points. Close to the impurity positions, the ef- 
fective interaction, comes mainly from Jh interaction, is 
ferromagnetic and becomes antiferromagnctic at a dis- 
tance comparable to ~ 2a. As expected, the amplitude 
of oscillations decrease with increasing distance between 



the local moment impurities. Moreover, the symmetries 
of the band structure which is a three-fold symmetric 
reflects in the susceptibility and thus in the magnetic 
interactions. The three-fold symmetric is clear in this 
figure. 

For a situation where the Fermi energy has an intersec- 
tion to the valance bands (hole doped case), the RKKY 
interaction terms exhibiting oscillatory behavior struc- 
tures however, the magnitude of the interactions are dif- 
ferent from those results obtained for the case of the elec- 
tron dope system. An intriguing manner of the RKKY 
interaction is its dependence of the Fermi energy values. 
The Fermi energy dependence of the RKKY interaction 
is shown in Fig. [4] for which the Fermi energy always in- 
tersects the conduction band and for a certain value of 
distance between two impurities, R = (4y3a, 0). The 
interaction amplitudes increase by increasing the Fermi 
energy and the period of the oscillations remain constant. 
It is expected that the values of J dm and J/, in the 
short-range, be small since the contribution of the spin 
splitting in the conduction band is negligible. Owing to 
a phase difference between Jj's, for certain values of the 
Fermi energy, there are situations in which J dm > Jh- 
Interestingly enough, for a certain values of the Fermi en- 
ergy and the distance between two magnetic impurities, 
J dm > Jh and leads to a strong rotation between spins 
of impurities. Accordingly, by designing an arbitrary ar- 
tificial lattice with magnetic impurities (e.g. triangular 
or square with arbitrary lattice constant), by properly 
adjusting the distance R using the scanning tunneling 
microscope techniques, we will be able to control each 
coupling term in the RKKY interaction and therefore 
the monolayer M0S2 would be a quite good platform to 
exhibit different spin lattices and interactions. 

On the other hand, for an undoped M0S2 system, the 
response of electron decreases exponentially and there- 
fore the RKKY interaction terms fall off rapidly. This 
phenomena is the same behavior which occurs in the 
topological insulators when the Fermi energy lies in a 
gapp2l. Fig. [5] shows the RKKY interaction terms for 
an undoped monolayer M0S2. Those interactions decay 
fast which denote that the interactions are rather short- 
rangecP^ and tend to zero quickly. Moreover, Jh and 
Ji are mostly positive which shows that the coupling be- 
tween the moments remains antiferromagnetic-like for all 
Rs. 

In the monolayer M0S2, there are two valance bands 
which are separated out by the spin-orbit interaction and 
furthermore, each band has a different s z value ( see 
Fig. 1). The top valance band starts at s = —A/2 + A 
whereas the another band starts at e = —A/2 — A. 
Fig. [6] shows the behavior of the RKKY interaction 
terms with respect to the Fermi energy in the region 
where the valance band is occupied. For the case that 
£f < —A/2 — A, an important result will be occurred. 
First of all, the RKKY interaction terms exhibit oscil- 
latory structures and their amplitudes suppress by in- 
creasing the Fermi energy. This behavior occurs till the 
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FIG. 2: (Color online) Long range behavior of RKKY inter- 
action terms, J; for i = H, DM and /, times to 7? 2 scaled by 
J c 2 /(27r) 3 a 2 as a function of the impurities distance along x 
direction for the Fermi energy which intersects the conduction 
band at ef = A/2 + A. The behavior of all J; interactions are 
determined by a superposition of two modified Bessel func- 
tions (see the text) with two different periods determined by 
arguments in j3±. Moreover the RKKY interaction terms in 
doped M0S2 decay as R~ 2 at long-range distance. 
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FIG. 3: (Color online) Short range contour plot of RKKY in- 
teraction terms, Ji for i = H, DM or I scaled by J 2 /(27r) 3 a 2 
as a function of the impurities distance in (x, y) plane scaled 
by a for the Fermi energy £f = A/2 + A. One of the spins is 
fixed at origin and another can be located on lattice points. 
Solid lines denote boundaries with positive values of the in- 
teraction term while dashed lines refer to its negative values. 
The results show clearly the three-fold symmetric. 
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FIG. 4: (Color online) RKKY interaction terms, Ji for i = 
H,DM or I scaled by J 2 /(27r) 3 a 2 as a function of electron 
energy doping in units of eV at R = (4\/3a, 0). Since there is 
a phase difference between JVs, for certain value of the Fermi 
energies, there are situations in which J dm > Jh- 



interaction terms decay rapidly. The reason for such de- 
caying comes from the fact that the one spin channel, for 
the such Fermi energy value, is blocked and therefore, the 
magnetic spin can not be coupled to the host electrons. 
In other words, in this case, the spin susceptibility of the 
system tends to zero since the spin relaxation time of 
the monolayer M 0S2 is quite long^ when the intervalley 
scattering is ignored. 
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FIG. 5: (Color online) Short-range behavior of RKKY in- 
teraction terms, Ji for i — H, DM or / scaled by J 2 /(27r) 3 a 2 
as a function of the impurities distance along x direction for 
undoped monolayer M0S2 system where ef = 0. Ji coupling 
interactions tend to zero quickly and moreover, Jh and Ji are 
mostly negative and the interaction is an antiferromagnetic. 

Fermi energy reaches to the top valance band illustrated 
by a shadowed layer in Fig. [6] and afterwards the RKKY 



FIG. 6: (Color online) RKKY interaction terms, Ji for 
i — H,DM or / scaled by J 2 /(27r) 3 a 2 as a function of the 
hole energy dope where e F < -A/2 — A for R = (4 v / 3a,0). 
The RKKY interaction terms oscillate and their amplitudes 
suppress by increasing the Fermi energy till it reaches to the 
top valance band illustrated by a shadowed layer. 



A. Twisted angle: A classical picture 

As discussed before, the DM term causes a twisted 
interaction between two magnetic adatoms. In order to 
calculate the equilibrium angel of the interaction, we sim- 
plify the RKKY interaction Hamiltonian in the spherical 
coordinate as: 

ilRKKY = hh(JH + Ji) cos 6>i cos 6 2 
+I1I2 sin^! sin 8 2 (J H cos <j) + J dm sin0) (12) 

where $1(2) denotes the angle between spin along z 
direction and <fi = cj>2 — 4>i- The equilibrium con- 
ditions, &H.rkky /99i = 0,dTiRKKY/d8 2 — and 
9Hrkky/9<P = 0, lead to the following equations: 
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{Jh + Ji) cos 82 sin B\ + sin 62 cos 9% {Jh cos < 
-{Jh + J/) cos 0i sin # 2 + sin 9\ cos # 2 (-?H cos ( 



sin 61 sin 62^— J h sine 



- Jdm sin( 

- J DM sin ( 

J DM COS ( 



= 
= 
= 



(13) 



There are solutions for Eq. ( 13 ) at the boundary of 
the square (9± — 8%) plane which are conveniently given 
by 7r(l,0), 7r(0,l), tt(0,0) and tt(1,1). For those so- 
lutions at the boundaries, possible configurations of 
the two magnetic impurities are only ferromagnetic- or 
antifcrromagnetic-like configurations. Moreover, there is 
just one extremum point at 9\ = 9 2 = 9 = n/2 be- 
tween 9 = and tt. The states of the system are de- 
scribed by a maximum, minimum or saddle behaviors 
of the extremum. An intriguing case can be occurred 
for the situation that the extremum is a local minimum. 
In this case, the equilibrium condition leads to the case 
that the moments lie in the sample's plane with an angle 
4> = tan _1 ( Jdm/Jh) between two magnetic moments. 

In order to describe the extremum behavior at point 
9 = it/2, we do need to evaluate the positive definite 
condition of the Hessian matrix of Eq. (12) throughout 
a neighbor of 9 = 7r/2, which is the second derivative of 
the RKKY Hamiltonian with respect to free variables^, 
and it yields as 



Di 
D 2 



Jdm ~ Ji(2Jh 4 



T 2 

J DM 



Ji) 



For a case that D\ > and D2 > 0, the critical point 
9 = tt/2 is a local minimum of the system of which we 
would like to explore the system. However, for the sit- 
uation in which Di < and D2 > the critical point 
is a local maximum of the system. A saddle point, on 
the other hand, is described by condition that £>i > 
and D 2 < or D\ < and D 2 < 0. For two latest 
cases, the equilibrium state of the system is obeyed by 
the boundary solutions which given by ferromagnetic- or 
antiferromagnetic-likc configurations. 

Figure [7] shows the sign of Z?i and D 2 as a function 
of the Fermi energy in the conduction band when the 
distance vector between impurities is R = (4-\/3a, 0). A 
saddle point, minimum and maximum regions are taken 
place by increasing the Fermi energy from the bottom 
of conduction, respectively. The inset figure shows the 
angle <f> between moments for the middle region in which 
both D\ and D2 are positive and 9 = tt/2. 



IV. SUMMARY AND CONCLUSIONS 

We study the influence of spin-orbit coupling on the 
RKKY interaction in a monolayer M0S2 system. We use 
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FIG. 7: (Color online) Sign of Di and D2 as a function of 
the Fermi energy in the conduction band for R = (4\/3a, 0). 
Inset figure shows the angle <j> = tan - {Jdm/Jh) between 
two magnetic moments when D\ and D2 are both positive. 



the Green's function method in the continuum model. 
We show that the rotation of the itinerant spin causes a 
twisted interaction between two magnetic adatoms which 
consists of different RKKY coupling terms, the Heisen- 
berg, Dzyaloshinsky-Moriya and Ising interactions. We 
explore different scenarios of two magnetic moments in 
terms of the different RKKY coupling values. We also 
find that the interaction terms are very sensitive to the 
Fermi energy values and change dramatically from doped 
to undoped systems. For undoped system, we find a 
short-range interaction between the moments with most 
probably an antiferromagnetic-like structure. The de- 
pendence of the interaction on the distance R between 
two local magnetic moments is found to be R~ 2 in a 
doped monolayer M0S2 sheet. We obtain that for a di- 
lute electron doping, the Dzyaloshinsky-Moriya interac- 
tion can be larger than the Heisenberg interaction for 
certain value of R. A beating pattern of oscillations of 
the RKKY interaction which occurs for the doped system 
is discussed. 

We conclude that the monolayer M0S2 is a quite good 
platform to explore different types of spin Hamiltonian. 
The interaction terms in the single layer M0S2 are quite 
similar to those terms appearing in the topological in- 
sulator systerrPH, however the interaction terms here are 
direction dependent where in the topological insulators 



the interactions are only distance dependent. 
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